Surface polaritons in two-dimensional left-handed photonic crystals 
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Using an extended plane-wave-based transfer-matrix method, the photonic band structures and 
the corresponding transmission spectrum of a two-dimensional left-handed photonic crystal are 
^jQ \ calculated. Comparisons between the periodic structure with a single left-handed cylindric rod are 

f"^ , made, and many interesting similarities are found. It is shown that, due to the localized surface 

f"^ ■ polaritons presented by an isolated left-handed rod, there exist many exciting physical phenomena 

£Nj ' in high-dimensional left-handed photonic crystals. As direct results of coupling of the localized 

, | , surface polaritons of neighboring left-handed rod, a lot of almost dispersionless bands, anti-crossing 

behavior, and a zero n gap are exhibited in the left-handed periodic structure. Moreover, in a 
certain frequency region, except distorted by a lot of anti-crossing behavior, there exists a continual 
dispersion relation, which can be explained by the long-wavelength approximation. It is also pointed 
■ out that high-dimensional left-handed photonic crystals can be used to design narrow-band filter. 
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Starting from the pioneering papers of Yablonovich [lj and John in 1987, the past decade has witnessed a great 
deal of interest in photonic crystals (PCs), also named as photonic band gap (PBG) materials. These have dielectric 
constituents which are periodically arranged in space, either in one, two, or three dimensions, with important potential 
applications in the optical, microwave and infrared fields . One of their important properties is to mould and control 
the flow and distribution of the light at its most microscopic level. The synergetic interplay in these PCs between the 
microcavity resonances of composite particles and the Bragg scattering resonances of the dielectric arrays lead to the 
"t- < formation of a PBG, i.e., a range of frequencies for which no propagating electromagnetic (EM) modes are allowed. 
O ,■ Due to the presence of the PBG in the dispersion relation of the EM field, the photonic density of states in PCs is 
'— 1 \ suppressed over a certain frequency window. These features open the possibility for many important technological 
applications including lossless PC waveguides, low-threshold PC lasers, and high-Q PC nanocavities 0. 

Most of the early research efforts were concentrated on PBG materials consisting of positive and frequency- 
. independent dielectrics. In the past few years^due to the new developments in structured EM materials, left-handed 
(LH) meta- materials have been fabricated 0, • This new kind of meta-materials were first speculated by Veselago 
in 1967, which characterized by simultaneous negative dielectric permittivity and negative magnetic permeability and 
therefore a negative refractive index in some frequency ranges 6]. Many anomalous effects can exist in such meta- 
materials, such as a reversed Doppler shift, reversed Cerenkov radiation, negative radiation pressure, and inverse 
Snell-Descarteds law. One of the most exciting applications of LH meta-materials is the possibility of a perfect lens 
whose image resolution is not subject to the limitations of the traditional diffraction limit of a wavelength |7J,|8(- 
Recently, a periodic arrangement of LH meta-materials, namely LH photonic crystals, have become the object of 
• i-H ' intense experimental and theoretical interest 0, 0, El El 0] ■ Like polaritonic PCs whose fi is homogeneous while e 
is frequency-dependent |l3j . the substitution of LH meta-materials into PCs is expected to introduce a whole range 
of exciting physical phenomena. For instance, a multilayer structure with stacking alternating layers of ordinary 
O (positive-n) and negative- n materials can result in an omnidirectional zero-n gap |9j, where the defect modes are 
L* ■ weakly dependent on the incident angles Jl4|. In general, the electric permittivity and magnetic permeability in LH 
meta-materials are frequency-dependent [8j. In a split ring resonator used to create negative magnetic materials, the 
magnetic permeability is normally expressed as n(u>) — 1 — ui^jioiui + ij), where u> m is the effective magnetic plasma 
frequency, u is the frequency of the light and 7 represents the light absorption Q . Due to the dispersive properties of 
LH meta-materials, it brings a lot of difficulties to simulate high-dimensional LH PCs with the conventional method, 
such as plane- wave expansion method 3] and finite-difference time-domain method 15] . 

Here, using an extended plane-wave-based transfer-matrix method, whose details are presented in the appendix, the 
optical properties of high-dimensional LH PCs are studied. The fundamental physical background of our method is: 
(1) Divide the whole PC structure into thin slices along the z-axis direction, within each slice e and fi are approximated 
as ^-independent (they can be xy-dependent); (2) It is assumed that each slice is surrounded by two thin air films 
with zero thickness; (3) Expression of the EM field in the air films by Bragg waves; (4) Expression of the EM field in 
each slice by the eigenmodes in the slice; (5) Use of the boundary conditions between EM fields in the slice and two 
neighboring air films to obtain the transfer matrix or the scattering matrix; (6) Use of the obtained matrix to calculate 
the transmission spectrum of the PC. The photonic band structure of the PC can thus be obtained by using Bloch's 
theorem. Since the independent variable in these calculations is frequency rather than wave vector, it can therefore 
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FIG. 1: The dependence on frequency uia/2nc of the parameters of the left-handed cylindric rods: fx(uj) (dotted), e(cu) (dashed), 
|n(w)| (dash dotted). 



effectively simulate LH PCs, even their magnetic permeability and dielectric permittivity are all frequency-dependent 
and negative. 

Based on this approach, the photonic band structures and the corresponding spectrum of two-dimensional (2D) LH 
PCs are calculated. The 2D PC considered here is a square lattice of LH cylindric rods embedded in air background. 
For the sake of easy comparison, the electric permittivity e(u>) and the magnetic permeability fi(uj) of the rods are 
adopted from Ref.jlg- The free electron expression of e(uS) is assumed, 



e( W ) = !--§, (1) 



to 



where lo p is the frequency of bulk longitudinal electron excitations, and 

, \ Flo 2 

MM = 1 - 2- 2 
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The parameters assumed are: F = 0.4, uj p = 0.552(27rc/a), uj = 0.276(27rc/a), and the radius of the rod r = 0.3a, 
where a is the lattice constant, c is the speed of light in vacuum. Note that the parameters loq and ui p used here are 
1/20 of that used in Ref.0. 

In Fig. (1) we plot the dependence on frequency Lua/2irc of the parameters of the LH cylindric rods: /Lt(w), e(u>), 
\n(co)\. It is clearly shown that: for frequency u> < u>q, e(w) < and n(u>) > 0, which means in this frequency region 
the meta-material is single-negative; for loq < u) < 0.35635(27rc/a), both and e(u>) are negative, that is, the 

meta-material is left-handed; for 0.35635(27rc/a) < to < u> p , e(w) < and fJi(u>) > 0, which means the meta-material 
is single-negative again; for uj > u) p , both n(w) and t{ui) are positive, and the meta-material is right-handed. Notice 
that 0.35635 is also 1/20 of the w& used in Ref. [l6T| . 

The photonic band structures for the transverse magnetic (TM) mode along the Y — X direction are calculated and 
shown in Fig. (2), and the corresponding transmission spectrum along the same direction for a 16-layer structure are 
also calculated and shown in Fig. (3a). Here, in frequency region u> < ojq and uj > 0.35 6 35(27rc/a) where the meta- 
materials are not LH, 101 plane waves are used, the cylindric rods are divided equally into 20 parts, and the increment 
frequency is 1.0 x 10 -5 (27rc/et); while, for the left-handed region, that is, u>o < uj < 0.35 6 35(27rc/a), 161 plane waves 
are used, the cylindric rods are divided equally into 50 parts, and the increment frequency is 1.0 x 10 _6 (27rc/a). It can 
be expected that, by decreasing the increment frequency, which leads to more computational burden, more accurate 
results will be obtained. 
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FIG. 2: The photonic band structures of the 2D PC. The inset shows the unbroken band structures. 




FIG. 3: Transmission spectrum of the 2D PC. (a) The frequency range is [0.2(2nc/a), 0.7(2nc/a)]. (b) The frequency range is 
[0.25(27rc/a), 0.37(2vrc/a)]. (3)The frequency range is [0.327(2?rc/a), 0.357(27rc/a)]. 



Comparing with the scattering resonances of a single LH cylindric rod shown in Fig. (2b) in Ref. ly|, there are two 
significant similarities between our periodic structure and the single rod :(1) Many anti-crossing behavior exist and 
only exist in the band structures where the meta-materials are left-handed, which lead to the formation of very narrow 
band gaps and the appearance of sharp resonant peaks in the transmission spectrum. Similar results were obtained 
from the scattering resonances of the single LH rod |16|. where resonances only exist in the frequency region where 
the cylindric rod is LH. The physical mechanism of the anti-crossing behavior in our periodic structure is same to 
that of the scattering resonances of the single rod, that is, the appearance of infinite surface polaritons [l6L Il7j: (2) 
A narrow band gap exists around 0.6(27rc/a), corresponding to the valley in the transmission spectrum. Notice that 
similar dip was also founded around the same place in the scattering resonances of the single rod 16]. It is a direct 
result of the fact that the effective refractive index of the periodic structure around 0.6(27rc/a), which is roughly given 
by 



n(u) = [ttt 2 ^e(uS) h(uj) + (a 2 - irr 2 ) x 1.0]/ 'a 2 , 



(3) 
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FIG. 4: The dispersion curves of our 2D LH PC (dotted) and the corresponding effective structure (solid line). 



is only 0.8. This value is very close to the refractive index of vacuum, therefore, the transmission rates around this 
narrow band gap is not zero, but still about 37%, as shown in Fig. (3a). Based on these two similarities, we can 
conclude that for a high-dimensional LH PC, just like a polaritonic PC, the band structure is governed primarily by 
the surface polaritons of a single rod, whose frequencies are determined by the rod's geometry, where the periodicity 
of the crystal will only introduce slight dispersion |13| . 

In addition, four more significant features can be found from the band structures of our periodic structure: (1) 
A wide band gap, from to 0.2545(27rc/a), is shown in Fig. (2) and Fig. (3a), which is mainly due to the single- 
negative property of the meta-material, so EM waves decays strongly when propagating through this structure; (2) 
As shown in Fig. (2) and Fig. (3b), there exists a continual dispersion curve in the frequency region [0.2545(27rc/a), uq] 
where e(uj) < and > 0. It is not so surprise to find that similar dispersion curves were reported in 2D 

polaritonic PCs, due to the fact that the meta-materials in this frequency region have the same optical properties as 
the polaritonic materials [l3|: (3) A lot of near-dispersionless bands appear when the meta-materials are LH, which 
are direct results of the coupling between the localized surface polaritons of an isolated cylinder. The most striking 
feature of the near-dispersionless bands is that their group velocities are almost zero. (4) From 0.32728(27rc/a) to 
0.35635(27rc/a), where the meta-materials are LH, the dispersion curve seems continual, except distorted by a lot of 
anti-crossing, as clearly shown in Fig. (3c) and Fig. (4). The dispersion relation here has the same physical origin as 
that in [0. 35635 (2tvc/cl), oj p ], where the meta-material is single-negative with a negative e(u>) and a positive /J,(u>). A 
direct evidence can be found from the unbroken band structure shown in the inset of Fig. (2), where the dispersion 
curve is continual from 0.32728(27rc/a) to about 0.6(27rc/a). This phenomena can also be explained by the long- 
wavelength approximation, that is, when the wavelength of the propagating EM wave is much longer than the lattice 
constant of PC, the dispersion relationship becomes approximately linear, and PC is effectively homogeneous under 
this circumstance. Note that the effective refractive index of the LH meta-materials in this region is very close to 
zero, as shown in Fig. (lb). In order to explain it more clearly, we plot the dispersion relationship of the effectively 
homogeneous structure in Fig. (4). The effective refractive index of the periodic structure used here is coarsely given 
by 

n(u)) = [-nr 2 ^e{uj)n{uj) + (a 2 - nr 2 ) x 1.0]/a 2 . (4) 

The dispersion relationship can be obtained by using the relation that k = n(Lo)uj/c. When plot it, we shift the 
frequency upwardly 0.0258(27rc/a). It is clearly shown that these two curves have similar slopes, which means the 
speed of the energy propagation (given by cko/dk) is almost same in the effective structure and that in the periodic 
structure. Hence, we can conclude that the long-wavelength approximation is validated here, and the group effective 
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FIG. 5: Transmission spectrum of the 2D PC around uja/2'KC = 0.293125. 



refractive index can be calculated by using Eq.(4). Note that Eq.(4) is very different with Eq.(3): the index obtained 
from last equation is always positive, in sharp contrast, the index calculated by the first one can be zero, or even 
negative. 

When the average refractive index of a multilayer structure n = J Q ° n{x)dxja equals zero, where a is the lattice 
constant, a band gap, namely zero n gap, will appear in the band structure 9J. Our 2D sample has also such zero 
n gap around u)a/2nc ~ 0.293125, as shown in Fig. (5). When the refractive index of the cylindric rods around this 
frequency is about —2.5368, the average refractive index of the 2D PC, which is given by Eq.(4), is zero. Clearly, 
in our case, the zero-n gap is only one of many mini gaps induced by the anti-crossing behavior. Moreover, around 
0.29307275(27rc/a), a very sharp peaks with very narrow bandwidth is shown in Fig. (5), and the ratio between the 
peak frequency and the FWHM of the peak is as little as 6.8 x 10 -8 . This property can be used to design narrow-band 
filter. 

In conclusion, using an extended plane-wave-based transfer-matrix method, the photonic band structures and the 
corresponding transmission spectrum of a two-dimensional left-handed photonic crystal are calculated. Comparisons 
between the periodic structure with a single left-handed cylindric rod are made, and many interesting similarities are 
found. It is shown that, due to the localized surface polaritons presented by an isolated left-handed rod, there exist 
many exciting physical phenomena in high-dimensional left-handed photonic crystals. As direct results of coupling 
of the localized surface polaritons of neighboring left-handed rod, a lot of almost dispersionless bands, anti-crossing 
behavior, and a zero n gap are exhibited in the left-handed periodic structure. Moreover, in a certain frequency 
region, except distorted by a lot of anti-crossing behavior, there exists a continual dispersion relation, which can be 
explained by the long-wavelength approximation. It is also pointed out that high-dimensional left-handed photonic 
crystals can be used to design narrow-band filter. 
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fund of Chinese National Science Foundation (10234040), Key Fund of Shanghai Science and Technology Foundation 
(05DJ14003). 



I. APPENDIX: THE EXTENDED PLANE- WAVE EXPANSION TRANSFER-MATRIX METHOD 

Divide a general three-dimensional (3D) PC into thin slices along the z-axis direction, within each slice e and /i are 
approximated as z- independent. Therefore, a partitioned lamellar slice can be approximated as a 2D grating |l8j . We 
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denote the primitive lattice of the grating by two unit vectors ai and a2, the corresponding reciprocal lattice by bi 
and h>i. We let an incident plane wave propagate along the z-axis with a wave vector k = {k 0x , fc 0y , ko z ). The EM 
field at an arbitrary point r can be expressed as a superposition of Bragg waves: 



E(r) = ^¥ Hj {z)e i{ki i'* x+ki *-v v \ 



(5) 



H(r) ^^Hytzje'^" 1 ^" 9 ', 



(6) 



where the Bragg wave vector = (fcij, x , hj^) = (ko x ,k 0y ) + ibi + jb 2 , Eij(z) and Hy(z) are expansion coefficients 
of the EM field to be determined. 

In the following we present a detailed description on how the eigenmodes of the EM field within a 2D lamellar 
grating slice are obtained on the plane-wave basis, and how the transfer-matrix connecting fields at two sides of the 
grating slice is related to these eigenmodes. We start from Maxwell's equations 



V x E(r) = zfcoM(r)H(r), 

V x H(r) = -ifc e(r)E(r). 



(7) 



Here e(r) and /i(r) are periodic functions of the grating slice, they are homogeneous along the z axis within the slice. 
We can rewrite the above equation into six partial-differential equations for (E x , E y , E z ) and (H x ,H y ,H z ). E z and 
H z can be deleted from these six equations, since the physical parameters within each lamellar slice are independent 
of z. This leads to the following four coupled equations 
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—H x = — — 

dz x ik n dx 
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(10) 



?-H = —- 

dz y iko dy 



H {dx y 



dy 



+ ik eE x . 



(11) 



We can further write down the plane- wave expansion expressions for dielectric permittivity and magnetic permeability, 



e ( r ) = E Ci J e 



iGij r 



(12) 



Gii r 



(13) 



(14) 



M" 1 (r)=E^ 1 e i 



iGa-r 



(15) 



Here = ib\ + jb 2 . Substituting them into Eqs. (4-7) yields 



^, ; .. = ^E ( WU,,, - k mn , y H mn , x ) + ik £ /WW 

mn mn 



dz 
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(16) 
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(19) 



We here introduce column vectors £=(•••, Eij, x , Eij iV , ■ ■ ■) and H = (■ ■■ , Hij^ x , Hij jV , • • •) , so Eqs. (12-15) can be 
written into a concise matrix form 



OZ 

d 

OZ 



where the matrices T\ and T2 are defined as 



\kij,y e ij:mnkn 



—h - f^ 1 h 



(20) 



(21) 



rpij;mn _ j_ ( ^ij,x^ij- mn kmn,y k,i^ x [l-. mn k mn ^ x k\(.ij- ;mn 

^0 \">a€ij- t mn ~ kij,yt J 'ij:mn"' mn >y ' t ij,yt J 'ij:mn' tmn ^ 



(22) 



It is easy to see that replacing fiij ;mn or fi i j. mn by Sij- mn (thus corresponding to a position-independent /i), the 
above matrices reduce to the ones in Ref. [Tsj . 

From the two matrices of Eq. I|21I22|) . the eigenmodes in the slice, expressed by eigenvalues 1 and eigenvectors 
S a , are obtained by solving — P = —T\Ti. Here we only choose (3 with Im(/3) > 0, which means that the for- 
ward propagating wave is absorbed by the media, instead of being amplified, for the sake of simple mathematical 
treatment. The eigenmode corresponding to /?| is written as E[ = E^^z) + E^^z), where E^ ^z) = E^ e l,3i ^ z ~ Zi \ 

E az( z ) ~ E~ e~ l @ i ( z ~ Zi \ For the moment, Ef and E~ are unknown. Furthermore, we define column vectors 
E+ = [. . . , E'^ i (z) 1 . . .] T , and E~ = [. . . , E^ { (z), . . .] T . The electric field column vector E can now be expressed by 
the superposition of all the eigenmodes E — S a {E+ + E~ ) . The corresponding magnetic field column vector can be 
obtained through H = T~dE/dz = T a (E+ - E~). 

A matching of the boundary conditions of the EM field between the grating slice and two neighboring air films, 
yields the transfer matrix 



T = 



an a u 
0-21 0,22 



Jf3h Q 





an a\2 
a>2i 0,22 



(23) 



where h is the thickness of the slice, an = (S~ 1 Sq + T~ 1 Tq)/2, ai2 — (S~ 1 Sq — T~ 1 Tq)/2, and 021 = ai2, 022 = an. 
So and To are the counterparts of S a and T a in an air film. So is a unit matrix, and To is a block-diagonal matrix. 

Since e~ l ^ h may lead to numerical instability for backward propagating waves [l8| . we use the scattering matrix 
here, which can be derived in the following way. We rewrite the above equation as 



( an 


aiA 




V «21 


a 2 2 J 





an ai2 \ ( Sli-i 
a2\ a 2 2 l lor, 



where fi^ and Qi-i are the electric fields in the two air films. 

a u nf + a 12 fir = e m ann+_i + e^ h a l2 nr_^ 

ai2^+ + ail Slr = e^a^tx + e~ m ■o u fir_ 1 . 



(24) 

(25) 
(26) 
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In order to avoid e l/3h , Eq. (22) is re-formulated as 

e m a 12 Q+ + e m a n ilt = a 12 nf-i + Hi^i-v (27) 

We thus arrive at 



an -e^a 12 \ / fi+ \ /e^a u -a 12 \ /Oti 
-e^ai2 an M^V-J l-ou e'%]U: 



(28) 



It can be shown that the inverse matrix of the left coefficient matrix has a form similar to 

( Pl P A (29) 
\P2 Pi J 

so that 

a nPl - e^ h a 12 p 2 = I, (30) 



aup 2 - e m a 12P i = 0. (31) 

From the second equation, we obtain p 2 = a u &® h a\ 2 p\. Combining with the first equation, the scattering matrix is 
finally obtained 

S= fpih+P2t 2 P lt 2 + P 2tl\ /g 2 x 
\Pit 2 +p 2 h p\tx+p 2 t 2 )' 

where p\ = [an - e^ h ai 2 a^ 1 e lf3h 'a^] -1 , p 2 = a^ 1 e zph 'ai 2 [an - e lf3h a 12 a n 1 e^ h a 12 }~ 1 , t\ = e i,3,l aii, and t 2 = -a i2 . 

Notice that the input of our approach is the frequency of the incident EM field. The method can therefore 
effectively simulate dispersive PCs, even when their magnetic permeability and dielectric permittivity arc all frequency- 
dependent. Furthermore, it is important to note that our method can also simulate defective PCs [Ifi, l2fl |. such as 
PC-based waveguides and microcavities, by combining with the supercell method. We can also utilize the symmetry 
of the simulated structure to substantially ease the requirement of calculation capability as was already pointed out 
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